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Unstable quasi g-modes in rotating main-sequence stars 
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Abstract. This paper studies the oscillatory stability of uniformly rotating main-sequence stars of mass 3-8 Mq 
by solving the linearized non-adiabatic, non-radial oscillation equations with a forcing term and searching for 
resonant response to a complex forcing frequency. By using the traditional approximation the solution of the 
forced oscillation equations becomes separable, whereby the energy equation is made separable by approximation. 
It is found that the K-mechanism in rotating B-stars can destabilize not only gravity- or pressure modes, but 
also a branch of low frequency retrograde (in corotating frame) oscillations in between the retrograde g-modes 
and toroidal r-modes. These unstable quasi g-modes (or 'q-modes') hardly exhibit rotational confinement to the 
equatorial regions of the star, while the oscillations are always prograde in the observer's frame, all in contrast to 
g-modes. The unstable q-modes occur in a few narrow period bands (defined by their azimuthal index m), and 
seem to fit the oscillation spectra observed in SPB stars rather well. The unstable q-mode oscillation spectrum of 
a very rapidly rotating 8Mq star appears similar to that of the well studied Be-star y, Cen. The unstable q-modes 
thus seem far better in explaining the observed oscillation spectra in SPB-stars and Be-stars than normal g-modes. 
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1. Introduction 

With the introduction of new opacity data by Iglesias et 
al. (1992) and a later update by Iglesias & Rogers (1996) 
it was suggested that the origin of (3 Cephei pulsations in 
massive main-sequence stars and of the slowly pulsating 
B (SPB) stars (Waelkens ll991|l of lower mass is caused by 
the classical opacity valve mechanism (Eddington 1930) 
operating on non radial g-mode or p-mode oscillations. 
Early stability analyses with these new opacity results in 
non-rotating stars have been published by e.g. Kiriakidis 
et al. I|1992[) . Gautschy and Saio (1993) and Dziembowski 
et al. (tlimi) . 

Observations by Balona (1994) and Balona and Koen 
( 1994) showed an unexplainable lack of SPB stars in two 
open clusters which do contain many pulsating (3 Cephei 
stars. It was suggested that the lack of SPB stars in these 
clusters might be correlated with the rapid rotation of the 
cluster stars. Ushomirsky and Bildsten l)1998|) . in an at- 
tempt to explain Balona and Koen's observations in terms 
of rotational stabilization of g-modes, applied the tradi- 
tional approximation combined with a quasi-adiabatic sta- 
bility analysis to study the effect of rotation on the k in- 
stability mechanism. They did not find a decisive solution: 
rotation seemed indeed to stabilize some pulsation modes 
but to destabilize other g-modes. This work was followed 
up by Lee (|2001|) who used a fully non-adiabatic analysis 
with a truncated expansion of spherical harmonics to ap- 
proximate the low-frequency g-mode oscillations in rotat- 



ing stars. Lee found some rotational stabilization caused 
by the Coriolis coupling with higher degree spherical har- 
monic components, especially for retrograde g-modes in 
the inertial regime, but this seems not sufficient to explain 
the lack of SPB stars in clusters with rapidly rotating B- 
stars. Townsend fffllBl , also applying the traditional ap- 
proximation, tried to explain the lack of observed rapidly 
rotating SPB stars by the well known effect of the Coriolis 
force on g-mode pulsations: the confinement of the oscilla- 
tion to the equatorial regions < a/(2Sl s )) of the star, 
where /i = cost? and a is the oscillation frequency in the 
corotating frame (e.g. Savonije et al. 0*995). The lack of 
rapidly rotating SPB stars would then be caused by a se- 
lection effect, not by rotational stabilization of g-modes. 
However, more recent observations indicate that the two 
open clusters studied by Balona and Koen may after all 
contain a few SPB stars (Stankov et al. l2002f) . so that the 
problem may be non-existent. 



There is, however, another complication: the existence 
of unstable rotational 'q-modes' in rotating main-sequence 
stars whose instability has been hitherto neglected. These 
q-modes exhibit hardly any rotational confinement to- 
wards the equatorial region, in contrast to normal g- 
modes. In this paper we will study the stability of these 
quasi g-modes in comparison with normal (retrograde) g- 
modes and check how well unstable q-modes can explain 
the oscillation spectra of rotating B-stars. 
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2. Basic oscillation equations 

We consider uniformly rotating main sequence stars with 
mass M s between 3-8 M Q and denote their radius by R s . 
We wish to study the oscillatory stability of these uni- 
formly rotating B-stars by subjecting them to a perturbing 
periodic force to determine the resonant behaviour. We as- 
sume the star's angular velocity of rotation ri s to be much 
smaller than its break-up speed, i.e. (fl s /il c ) 2 <C 1, with 
fi 2 = GM s /i?g, so that effects of centrifugal distortion 
(oc SI 2 ) may be neglected to a first approximation. The 
Coriolis acceleration is proportional to Jl s and we consider 
its effect on the induced oscillatory motions in the star. 
We use nonrotating spherical coordinates (r, z9, ip), with 
the origin at the star's centre, whereby $ = corresponds 
to its rotation axis. 

Let us denote the displacement vector in the star by £ 
and perturbed Eulcrian quantities like pressure J", density 
p', temperature T' and energy flux F" with a prime. The 
linearized hydrodynamic equations governing the non- 
adiabatic response of the uniformly rotating star to the 
perturbing potential $t may then be written as 
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where k denotes the opacity of stellar material and S its 
specific entropy. These perturbation equations represent, 
respectively, conservation of momentum, conservation of 
mass and conservation of energy, while the last equation 
describes the radiative diffusion of the perturbed energy 
flux. For simplicity we adopt the Cowling (1941) approxi- 
mation, i.e. we neglect perturbations to the gravitational 
potential caused by the star's oscillatory distortion. For 
the higher radial order oscillation modes studied here, 
this approximation should be adequate. As usual, we also 
neglect perturbations of the nuclear energy sources (not 
important) and of convection (frozen convective flux ap- 
proximation). In main-sequence (MS)stars more massive 
than about 6 Mq in which the opacity Z-bump driving 
region is convectively unstable this latter approximation 
introduces some uncertainty with respect to the detected 
oscillatory instabilities. 

For the periodic forcing we apply a force equal to the 
(real part of) the gradient of 

$ T (r, 0, ip, t) = fr l PTifi) e^- m ^ (5) 

where a is the forcing frequency in the inertial frame, 
fx = cos$, -P ; m (At) is an associated Legendre polynomial 



and / is an arbitrary constant. Since we are interested in 
periodic solutions of the forced oscillations we assume the 
perturbations have the same time and ip dependence as 
the forcing (although with a certain phase lag due to ra- 
diative damping) . After factoring out the known time and 
(/^-dependence, Eqs. (^QJ form a 2D problem in the (r, d) 
meridional plane of the perturbed star. The 2D problem 
requires a large amount of computer memory and CPU 
time to solve, see e.g. Savonije et al. I|1995fl . However, it is 
possible to obtain solutions adequate for initial analyses 
of the stability of rotating stars by applying the so called 
'traditional approximation'. 

3. The traditional approximation 

For a rotating star the solutions of Eqs. are no longer 
separable into r-,i?- and ip- factors and cannot be described 
by a single spherical harmonic (l,m) due to the action of 
the Coriolis force on oscillating stellar matter. However, 
for a uniformly rotating star in the traditional approxima- 
tion separability is retained by neglecting the compo- 
nent of the angular velocity vector (e.g. Unno et al. 1*0)89). 
This is done because the radial motions are expected to be 
small in the stably stratified layers of the star, especially 
for low-frequency modes (Savonije et al. [1995). In this 
way the problem is reduced to two 1-dimensional prob- 
lems (in i? and r) with a prescribed azimuthal harmonic 
dependence e mip since we assume the unperturbed star to 
be axisymmetric. When only the radial component of Q s 
is retained, the i?- and ^-components of the equation of 
motion can be written: 



ct 2 ^ — 2ifi s CT cos i9£a = — 



a + 2iil s a cos 



1 dP' l<9$ n 



rp dd 



rp sin ■& 



:P' 



It) 



(6) 



(7) 



We have expressed the velocity perturbations in terms of 
the displacement vector by the relation v' = irr£, with 
a — a — m VL S the oscillation frequency in the corotating 
stellar frame (negative if the oscillation is retrograde, i.e. 
propagates in the direction counter to the stellar rotation). 
We will only consider positive m-values in this paper, ex- 
cept when explicitely stated otherwise (in Tables EEJ- 

3.1. Separation of variables 

A separation of variables can be performed by writing, see 
Papaloizou & Savonije l|1997f) : 



& = JL ^ (0) D n (r ) £ v = J2 Gn (0) D n (r) 

n=0 n=0 
oo 

P' = J2 X nWW n (r), (8) 
n=0 

with £ r , T 1 and p 1 having expansions of the same form 
as that for P'. Here we leave out the known factor e %mtp 
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for the azimuthal variation of the perturbed quantities, 
while T n , Gn and X n are functions of i? chosen to obey the 
relations 



— u T n — 2if2 s <7 COS 1? Q n 

and 

— <J 2 Gn + 2ifl s <7 cos d T n 



dXn 



sin$ 



We obtain an equation for X n ($) by imposing the con- 
straint 
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where A„ is a constant. Then X n must satisfy the second- 
order equation obtained from 
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(10) 



whereby is an eigenfuction with A„ the associated 

eigenvalue. The solution depends on the rotation param- 
eter x = 2Vl s /<j. For r2 s = the functions X n {"d) become 
the associated Legendre functions P™ +n (cos^) with cor- 
responding eigenvalues A„ = (to + n) (m + n + 1). Normal 
modes of the rotating star correspond to normal )modes 
of the non-rotating star with (to + n) (to + n + 1) replaced 
by any permissible A„. 

Different X n ($) are orthogonal on integration with re- 
spect to p = cos $ over the interval (—1, 1). If the perturb- 
ing potential is expanded in terms of the X n such that 



$tM) =X> n (r) 



we find from © and that 
D n (r) = — W n (r) + -*„(r). 



(11) 



(12) 



This equation is exactly the same as in the non-rotating 
case. The same is true for the equation of continuity, ex- 
cept that (m + n) (m + n + l) is replaced by A„. In this way 
the adiabatic stellar response will consist of a superposi- 
tion of responses appropriate to non-rotating stars with 
(m + n) (to + n + 1) replaced by A„ obtained from equa- 
tions H9I1()(I . However, we are interested in the oscillatory 
stability of rotating stars and for that we need to include 
non-adiabatic effects, i.e. we have to consider the energy 
equation . We can write the divergence term on its right 
hand side as: 



By adopting T'(r, §) = £ n T n (r) X n {<&) and substitution 
in the flux equation @ the radial part of the above di- 
vergence can be written in the required separated form 
T> r (r, •&) — H(r) X n ($). However, the two angular terms 
cannot be expressed in this way. For one can write, after 
applying equation 10} and equations l|9lll)[l . the angular 
part of the above divergence as: 



1 /dlogTA 1 T'(r) 



dr 



T(r) 



with 



3>n(0) 



2flx 2 Qn(p) , , .2 ,.„> 



+ ( 



X fj, — TO X 



-Xr 



which shows that the energy equation renders the sys- 
tem of equations inseparable, unless x — > 0. Nevertheless, 
we wish to apply the traditional approximation and take 
advantage of its simple treatment of rotational effects. 
Therefore we will approximate the energy equation by a 
separable equation by averaging the angular part of the 
divergence over $ and using the averaged value e n for y n ■ 

_ J y n X n dp 



J X n dp, 



so that 



1 fdlogT\ 1 T'(r) 
dr J T(r) 



X n {d) 



(13) 



For small values of the rotation-parameter x — > the exact 
result e„ — ► — A„ = —(m + n) (m + n+ 1) for non-rotating 
stars is approached. 

Expression (|13|) takes the radiative diffusion in the hor- 
izontal direction at least qualitatively into account. Below 
we will check our results by comparing with other linear 
stability calculations of rotating stars. 

3.2. The radial part of the oscillation equations 

Once we have solved the above angular eigenvalue problem 
for A„ we can factor out the common, now known, factor 
Xni^e)^*-" 1 ^ from all perturbed quantities. The lin- 
earized equations describing the non-adiabatic forced os- 
cillations, for the n-th component in expansion i|ll|) . then 
form a 1-dimensional (radial) problem. The remaining ra- 
dial part of the oscillation equations can be expressed as: 
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(15) 
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where T's represent Chandrasekhar's adiabatic coeffi- 
cients, 

^'diogry 1 ^ = dio g p r ^diogp 



dr 



dr 



the opacity derivatives are given by 
d log k\ ( d log k 



Kt 



V 



dr 



d\ogTj p \d\ogp /T 
rj is a characteristic radiative diffusion length: 

4acT 3 dT 



-(r 3 - l)F/ (aP) , with F 



3np dr 



being the unperturbed (radial) radiative energy flux, e n is 
related to the horizontal radiative diffusion and is defined 
by (|13|l . All other constants have their usual meaning. 
Note that the radiative diffusion introduces a factor i, so 
that the radial parts of the perturbations are complex- 
valued. This expresses the induced phase-lags with respect 
to the external forcing caused by non-adiabatic effects due 
to radiative diffusion. 

The oscillation equations are complemented by the lin- 
earized equation of state 



EL - 
P ~ 

with 

A4 = 

Pa 



<91ogP 
dlogT 



d lOg Pa 

dr 



dlogP 
dlogp 



d log P 
d log p a 



^(18) 

Ma 



where p a is the mean atomic weight of the stellar gas. 

Finally, we prescribe the usual boundary conditions 
£ r = F' = at the stellar centre and require that the 
Lagrangian perturbations at the stellar surface obey: SP = 
and 5F/F = AST/T (Stefan's law). 

3.3. Stellar input models for the oscillation code 

We constructed the unperturbed stellar models for the 
main-sequence stars with a recent version (Pols et al. H995fl 
of the stellar evolution code developed by Eggleton l|1972[l . 
The models represent spherical main-sequence stars with 
masses 4— 8M and chemical composition given by various 
values of the central hydrogen abundance X c and Z = 0.02 
or Z = 0.03. These models were constructed by using the 
OPAL opacities (Iglesias & Rogers ESS- TableUJlists the 
effective temperatures of the various stellar input models 
characterized by their mass, core hydrogen mass fraction 
X c and metal content Z, used in our calculations. 

4. Solution method 

4.1. Determination of the eigenvalue X n (x) 

For a given stellar rotation frequency Q s , forcing frequency 
in corotating frame <7, azimuthal index m and angular or- 
der n, the eigenvalues A„ can be determined by numer- 
ically integrating the two first order differential Eqs. © 



Table 1. ( 10 log(L/L ), 10 logT cff ) of the input models; an 
asterisk above stellar mass indicates models with Z=0.03 
instead of Z=0.02 



Model 


X c = 0.6 


X c = 0.4 


X c = 0.2 


3M Q 


(1.932,4.069) 


(2.025,4.037) 




4M 


(2.411,4.153) 


(2.517,4.124) 


(2.597, 4.0790 


6M Q 




(3.173,4.238) 


(3.267, 4.199) 


6*M 




(3.112,4.210) 


(3.205,4.167) 


8*M© 




(3.562,4.288) 


(3.662, 4.247) 



and (|10|l . We have used a shooting method with fourth or- 
der Runge-Kutta integration with variable stepsize (Press 
1992) to obtain numerical solutions. It is convenient 



et al 

to rewrite equations I|9I10I) as 



dX n 
dp 

dQn 

dp 
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1 
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1-p 2 



— X n I X n 



mxp 



1 



■li* 



Qn 



(19) 



(20) 



where the factor x — 2il s /a expresses the importance of 
stellar rotation. 

To enable integration away from p = 1 it is convenient 
to write X n (p) — (1 — p 2 )^Y n (p) and expand Y n (p) in a 
power series which can be substituted in Eqs. (|19|) and l|20|l 
to determine the coefficients. Q n can then be expressed in 
terms of Y n as 



Qn 



(i-M 2 y 



i 



x 2 p 2 



dY 

mp(l - x) Y n - (1 - p 2 ) -j^ 



For the even n (including n—0) solutions the boundary 
conditions are X n = at p = 1 and Q n = at p = 0, 
while the odd solutions are defined by X n = at p = 0. 
We can now integrate equations (fTl^l - (|271)l from p = 1 — 5 
(with 5 = 10~ 4 ) to p = with an estimated value for 
the eigenvalue A„. We iterate by adjusting A„ until the 
integrated value for either Q n (even solutions), or X n (odd 
solutions) is sufficiently close to zero for p = 0. 

4.2. Eigenvalues \ nm and mode classification (n,m) 

Eigenvalues X nm of equations I|9I10() (or its equivalent) 
have recently been calculated by Bildsten et al. ©96), 
Papaloizou & Savonije l|1997[) . Lee & Saio l|1997|l and 
Savonije & Witte ( 2002). It appears that in the traditional 
approximation the angular eigenvalue X nm can take both 
positive and negative values, including 0. Here we will fol- 
low Lee & Saio (1997, see their Fig. 1) and allow n to 
assume positive and negative values (or 0) to discriminate 
between these cases. The positive values n = 0,1,2,3,... 
correspond to gravity (g)-modes (or p-modes) for which 
A > for any x-value, whereby even values of n (or 0) 
correspond to modes with even symmetry about the stel- 
lar equator and odd values to modes with odd symmetry. 
For these modes n gives the number of nodes in X n (p) 
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Fig. 1. The eigenvalue A„ versus a/(2£l s ) for m=2 modes. 
From top to bottom right: full line: g-modes with n = 0; 
dashed line: q- modes with n = — 1; dotted line: q- modes 
with ri = —2. The n = — 1 branch of quasi g-modes lies 
between the normal g-modes and the quasi-toroidal modes 
(with A_2 remaining small as \l/x\ —> 0). The eigenfunc- 
tion X n {"&) corresponding to n = and n = —2 has no 
nodes, that of n = — 1 has a node at the equator. 

over the interval — 1 < fi < 1, except when — 1 < x < 
(where x" 1 — a /2f2 s ) when there appears an extra set of 
nodes. This is related to the existence of quasi g-modes 
with n = — 1, see below. 

For n < and x < the eigenvalues are positive for 
— » 0, and negative for large values of |x _1 |. In this 
case the eigenvalue A nm changes sign for x-value: 



(m+ H) (to + \n+ 1|) 

where to > 0. The oscillation modes corresponding to 
these roots the purely toroidal r-modes (Papaloizou 

& Pringle 1978). The n < branches thus have positive, 
although small, eigenvalues A„ m for |x _1 | < \x~ 1 \ and 
these correspond to quasi-toroidal q-modes (Savonije & 
Witte, 2002). Fig. [I] shows the (A > 0) q-mode branches 
with n — — 1 and n = —2 for m—2 together with the 
normal g-mode n—0 branch. The special case n = — 1 
is interesting as this corresponds to a branch of q-modes 
where A does not remain small when |x _1 | decreases, see 
Fig-E For | 1 1 < la;,: 1 ! these n = — 1 modes have values 
of A nm that remain an order of magnitude smaller than 
those of the lowest order (n — 0) g-modes. This means 
that for these special 'quasi g-modes' the confinement to 
the equatorial regions of the rotating star at lower fre- 
quencies a is almost absent. We will see that for sufficient 
rotation speeds n = — 1 quasi g-mode oscillations can be 
destabilized by the opacity bump due to the heavy element 
ionisation zone, like normal g-modes. 

The n < modes with X nm < (for both positive 
and negative x-values) correspond to rotationally stabi- 
lized g - -modes or oscillatory convection (and exists only) 
in the regions where the square of the Brunt- Vaisala fre- 



quency AT 2 = g A/Ti < 0. In this paper only modes with 
Anm > are explored. 

4.3. Solution of radial oscillation equations 

For a given stellar rotation speed fi s and stellar forcing 
frequency a (in fact for a given x^ 1 = a/2il s ) we can for 
each possible combination of (n, to) determine the eigen- 
value A n and corresponding eigenfunction X n (d) (and aux- 
illiary function Q n from which e„ is calculated) , substitute 
A n and e n in the radial oscillation equations (|15I16I) and 
solve the set of radial equations I|14ll7f) . The solution of 
these differential equations is found by transforming them 
into algebraic equations by means of finite differences on 
a staggered spatial mesh and applying matrix inversion 
similar to standard Hcnycy schemes for stellar evolution 
(e.g. Savonije & Papaloizou 1983). Hereby the forcing po- 
tential ^t(1,iti) is chosen to have the correct symmetry, 
i.e. for given (n, m), we adopt / accordingly: for g-modes 
I = m + n and for q-modes I = m+ | — n — 2\. By scanning 
through a range of (initially) real forcing frequencies a and 
solving the radial oscillation equations we find the possible 
oscillation modes (with a different number of radial nodes 
k) when the response becomes resonant for certain values 
afc(f2 s , n, to). At resonance (that is what we are interested 
in) the stellar response becomes that of the free oscillation 
mode (fc,n, to) and appears virtually independent (apart 
from a fixed overall scaling factor) of the value adopted 
for / in the forcing potential, as long as the equatorial 
symmetry of the forcing is consistent. 

4.4. Search for unstable oscillation modes 

The resonances with free oscillation modes are found by 
varying the real valued a and maximizing the absolute 
value of the tidal torque integral (interpreting the forcing 
potential as that of a spherical harmonic component of 
an external point mass companion as in Savonije & Witte 
2002): 

T nm (a, n s ) = it Cnm / / Im [p' nm (r)} r l+2 dr (22) 
Jo 

where 7r results after the integration over ip, f is the con- 
stant in the tidal potential defined by Eq. (JSJ), Im stands 
for imaginary part, and 

(The tidal torque follows by multiplying T nm with m). 
This procedure shows in fact immediately whether the 
resonance corresponds to a stable oscillation mode (if T nm 
has same sign as a) or to an unstable one (when T nm has 
opposite sign as a) . This can be understood physically by 
considering a stable oscillation mode: when the forced os- 
cillation is prograde ( a > 0) the torque should lead to a 
spin up of the star and must be positive. For an unstable 
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mode the phase difference between forcing and star has 
the opposite sign. 

For a check of the stability properties of a mode 
and in order to estimate its growth or damping rate, we 
forced the star with a complex frequency a and searched 
with a numerical algorithm for that value of the imag- 
inary part Xm(cr) for which the integral over the star 
J = Jq ' Cr(r)C*i r ) becomes maximized. The search in- 
terval for Xm(cr) was chosen symmetric about the zero 
value. The resonant forcing of stable oscillation modes be- 
comes enhanced when the radiative damping is neutral- 
ized, so that J becomes maximal for a positive value of 
Zm(er) (identified as the 'damping rate'), in analogy with 
the standard theory for a ID forced linear harmonic os- 
cillator with damping. For an unstable mode with 'nega- 
tive damping', the enhanced resonant response is found for 
Tm(a) < (the absolute value of which is identified as the 
'growth rate'). Except for modes near the boundaries of 
the instability bands, the integral J of the squared oscilla- 
tion amplitude becomes many orders of magnitude larger 
when applying complex forcing frequencies to search for 
the maximum resonant response. 

Please note that in this paper we consider a to be 
the real part of the forcing frequency, unless its complex 
character is made explicit, like in Tm(o~), while a always 
refers to the real (part of the) oscillation frequency (in the 
corotating frame). 

4.5. K-mechanism 

The instabilities found in this paper are (checked to be) 
caused by the well known k valve-mechanism by which 
a fraction of the outflowing stellar thermal energy-flux is 
tapped by an increase of the radiative opacity during the 
compression phase of the oscillations and heat is released 
during the expansion phase. According to a semi-adiabatic 
analysis assuming constant radiative luminosity the opac- 
ity k has the correct absorbing behaviour in regions of the 
star where (e.g. Unno et al.|l989): 

s(- + f^r)>° < 23 > 

Regions where inequality (|2*3)) is fulfilled tend to drive the 
oscillation when a sufficient amount of thermal energy can 
be absorbed during the compression phase. For this the 
oscillation period should be comparable to the thermal 
timescale in the driving zone. Let us define a local 'adia- 
batic frequency' v a d = 27r/rj e afc and use the leaking time 
Ti ea k — (R s — r) 2 pn/c as a crude estimate for the ther- 
mal timescale r t h at radius r based on random walk of 
photons towards the stellar surface. When the surface is 
approached the leaking time becomes ever shorter and u a d 
rises steeply, see e.g. Fig. 0] This figure shows the various 
frequencies as a function of the temperature in a 4 Mq 
main-sequence star with Z=0.02. 

Note that from now on we normalize all frequencies on 
the star's critical rotation frequency il c — -J G M s /B? s also 



when this is not indicated. In the surface regions where the 
oscillation frequency becomes much smaller than the 'adi- 
abatic frequency' v a d the oscillations evidently are very 
non-adiabatic and no driving can occur when lll'-il) attains 
positive values in these regions. For instability the rise of 
v a d, where its curve in Fig. 0] intersects the actual (con- 
stant) a should occur well into the driving region where 
(|23|l attains a significant positive value, i.e. in the opacity 
bump region near T ~ 2 x 10 5 K, or rise steeply early in the 
damping region adjacent to it (so that the extra damping 
has no effect due to the already strong non-adiabaticity). 

It is known (Berthomieu et al. 119781 Savonije et al. 
1995} that the traditional approximation yields qualita- 
tively correct results for g-modes if the oscillation fre- 
quency a is small compared to the Brunt- Vaisala fre- 
quency (thus outside the convective core), so that the 
horizontal oscillatory motion dominates. This is also a rea- 
sonable approximation for the q-modes with weak radial 
motion studied here. Fig.0|shows that this requirement is 
amply satisfied for the 4M Q star. 

4.6. Comparison with other stability analyses 

We applied our calculation method to a non-rotating star 
and compared the results with those given in Gautschy 
and Saio Q1996[l and found an instability region for a 5M Q 
MS star very similar to the one shown in their Figure 2. We 
also compared our results for rotating stars with those of 
Lee H2001fl who used a truncated series expansion in spher- 
ical harmonics to describe the non-adiabatic non-radial os- 
cillations in a rotating star. Lee applied this scheme to a 
4M Q ZAMS star with X c =0.7 and Z=0.02. We will com- 
pare the results for rotation speeds S1 S =0. 1 (Table[2} and 
il s = 0.3 (Table EJ. Our 4M Q ZAMS stellar input model 
has lo log( J L/L )=2.366 and 10 logT c ff=4.170 and thus is 
slightly hotter than Lee's input model ( 10 logT c ff=4.164). 
The results presented in Tables 12131 show that the two 

Table 2. Comparison with Lee's (2001) result for the ra- 
dial orders (k m i n ,k max ) of unstable g-modes for a 4M 
ZAMS star (Z=0.02), rotating at angular velocity f2 s =0.1. 
The upper pair (k m i n ,k max ) corresponds to the results 
obtained here, the lower pair are those derived by Lee. 
For easy comparison we have given the (l,m) values with 
I — n + \m\. 



degree 


m = — 1 


m = -2 


m = —3 


m = 1 


m = 2 


m — 3 


/ = 1 


(10, 14) 
(9,15) 






(9,23) 
(9,23) 






I = 2 


(9, 19) 
(6,20) 


(9, 19) 
(8, 20) 




(9,23) 
(9,22) 


(9,22) 
(9,22) 




I = 3 


(10, 24) 
(9,24) 


(10,23) 
(9, 23) 


(10,23) 
(9, 23) 


(10,24) 
(9, 24) 


(10,24) 
(9, 24) 


(10,24) 
(9,24) 



methods yield very similar instability intervals for a 4M Q 
ZAMS star. Only for the highest rotation rate we find 
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Table 3. Comparison with Lee's (2001) result for the ra- 
dial orders of unstable g-modes for a 4M ZAMS star 
(Z=0.02), rotating at angular velocity f2 s =0.3. 



degree 


m — —1 


m = -2 


m — —3 


m = 1 


m = 2 


m = 3 


I = 1 


(10,13) 
(4, 15) 






(9,26) 
(11,24) 






I = 2 


(9,22) 
(7, 23) 


(9, 19) 
(6,20) 




(9,26) 
(7,22) 


(9,26) 
(8,22) 




I = 3 


(10,26) 
(7,21) 


(10,25) 
(9,24) 


(10,22) 
(8,23) 


(10,27) 
(8,20) 


(10,27) 
(9,20) 


(10, 26) 
(10, 24) 



slighly narrower instability domains for prograde g-modes 
on one hand and slighly more extended (to lower frequen- 
cies) intervals for retrograde g-modes on the other. This 
may be partly due to the slightly different input model. 

It makes practically no difference to these results 
whether we adopt e„=0 or e„ ^ 0, in the former case 
the instability interval is, for the case fi s =0.3, sometimes 
extended whereby k m in decreases by one. The great ben- 
efit of the present calculation method is that there is no 
truncation problem in the solution of the angular eigen- 
functions (although by neglecting the horizontal compo- 
nent of O s , which seems reasonable for the considered 
modes). The horizontal part of the energy flux perturba- 
tion is treated only approximately, see equation (|13fl . But, 
since the approach in this paper greatly simplifies the anal- 
ysis of the stability of oscillation modes in rotating stars 
and provides us with an at least qualitatively correct pic- 
ture of the linear stability properties of these modes, it 
seems a good initial approach. Townsend (2005) used a 
similar approach (but neglected the horizontal diffusion 
altogether by adopting e„ = 0) to study the stability of 
g-modes in rotating SPB-stars in more detail than is done 
here. Where the results overlap they seem qualitatively 
consistent with each other, although Townsend finds un- 
stable n=0 g-modes even in ZAMS stars somewhat less 
massive than 3M© while the n=0 g-modes appear all sta- 
ble in our 3M© models, see below. This seems related to 
the higher effective temperatures of the lower mass stellar 
models used by Townsend as compared to the models used 
here. 

5. Results 

5.1. Unstable n = — 1 q-modes in a 3M Q MS star 
(Z=0.02) 

Before exploring the stability of more massive MS stars 
let us search for the lower mass limit for which overstable 
q-modes can exist. We studied the stability of two MS 
models with masses 3M© and 2M Q . For the 2M Q MS star 
no q-mode instabilities (nor any unstable n—0 g-modes) 
could be found, so that the lower boundary of the q-mode 
instability region lies between 2Mq and 3M Q , coinciding 
by the way with the lower boundary for SPB stars. It can 
be seen in Fig. that for these lower stellar masses the 



Z-bump in the opacity near 10 logT ~ 5.3 is rather mod- 
est and because it lies quite deep inside the star, where 
the local thermal timescale Tth is relatively long, only 
very low frequency oscillation modes can be sufficiently 
non-adiabatic in this potential driving zone to give rise to 
overstability. It appears that the very low frequency quasi 
g-modes (or n=-l 'q-modes') can indeed be excited in the 
3M star, see Table 0] for a list of unstable modes. No 
unstable g-modes could be detected in a 3M Q MS star: 
the required low frequencies can only be reached for high 
radial order modes for which internal damping connected 
with short wavelengths is too strong. 
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Fig. 2. MS star of 3 M© with A c =0.6 and Z=0.02. 
Full line: the modified Brunt- Vaisala frequency 
t'BV = sign(A/" 2 ) y/ 1 N 2 1 , where N is the Brunt- Vaisala 
frequency; dashed line: the 'adiabatic frequency' 
v a d = ^Tt/Tieak', dot-dashed line: again is a( i, but now 
for X c — 0.4; dotted line: the value of the radial deriva- 
tive in expression l|23|l in arbitrary units, all versus 
temperature. 

Fig. El shows the oscillation frequencies in the corotat- 
ing frame: ak/(2Q s ) of n=-l q-modes for several rotation 
rates of the X c = 0.6 and X c =0A model. The (absolute 
value) of the frequencies in the corotating frame |<7fc| are 
for m > 1 significantly lower than those of unstable normal 
g-modes. For given f2 s the mode frequencies |<7fc| are seen 
to decrease with increasing m- value, while for given m the 
frequencies \ok \ increase with rotation speed Q s . That cr^ 
increases with m follows from equation l|21() which tells us 
that the roots \x~ 1 \ from where \ nrn > for n = — 1 
modes shifts towards (respectively x^ 1 = —1/2,-1/3 
and —1/4 for m—1, 2 and 3), so that the lowest radial or- 
der mode frequency |cto| for a given m is larger than that 
for m + 1 and consequently [tTfc(m)| > \dk(jn + 1) for not 
too large k. 

For m > 1 the mode frequencies of adjacent radial 
orders k are very densely packed, so that the oscilla- 
tion periods of all q-modes in an instability interval 
(kmin, k max ) can only be distinguished if the observations 



8 



G.J. Savonije: Unstable quasi g-modes in rotating main-sequence stars 



1 4- 

- 


4- 


i 

+ 


1 

-t- 


i 

+ 


1 1 

- 


- x x 
■ 


X 
■ 


X 

■ 


X 
■ 


X 


* X 




i 


G 


* 
D 




* * 


i i 


i 


I 

1 


B 

i 


i 

1 


e 

i i 



13 14 15 16 17 18 19 20 21 22 

k 



Fig. 3. Unstable n=-l q-modes: (7fc/(2fi s ) versus the num- 
ber of radial nodes k. £l s =0.3: plusses: m=l, crosses: to=2, 
asterisks: m=3 and dotted squares: to=4. Then for f2 s =0.4: 
filled squares: m=2; dotted circles: to=3 and filled circles: 
m=4. Note that for given f2 s the n—-l q-modc frequencies 
\&k\ decrease with increasing to. Model: 3 M©, X c = 0.6, 
Z=0.02. 

have good time coverage, see Table 0J This property of 
q-modes tends to give a much simpler observable oscilla- 
tion spectrum in rotating stars than g-modes. Note also 
that although in the corotating frame the q-modes are 
retrograde modes, in the observer's frame they are always 
prograde. 

The location and extent of the instability intervals are 
very sensitive to the thermal timescale Tth in the opacity Z- 
bump region. In the more evolved MS models, with cooler 
surface layers, this region lies deeper inside the star so that 
T t h is longer and the steep rise of the 'adiabatic frequency' 
curve takes place a bit further outwards in the star, at 
lower T, see dot-dashed curve in Fig. |2l Instabilities in 
the less evolved X c =0.6 model with the shorter r t h re- 
quire somewhat higher oscillation frequencies for instabil- 
ity and can thus occur at higher rotation rates than in 
the more evolved X c =0.4 model, see Table 0| For rota- 
tion speed f2 s =0.3 instabilities can only be found in the 
A c =0.4 model for to > 3 (i.e. for the lowest mode frequen- 
cies) , while for still higher rotation rates (and thus higher 
mode frequencies) no q-mode instabilities could be found. 
We see in Table 0] that the modes with higher to values 
(i.e. with lower frequencies |crfc|) become stable when Q s 
decreases as the oscillations become too slow and hence 
too non-adiabatic. 

5.2. Unstable modes in a 4M & evolved MS-star 

5.2.1. Unstable retrograde n=0 g-modes 

Let us now turn to 4 M© MS star models of decreasing 
central hydrogen abundance X c . Fig. 01 shows the (mod- 
ified) Brunt- Vaisala, frequency, the 'adiabatic frequency' 



Table 4. Unstable n = — 1 q-modes for a 3M© star with 
Z= 0.02 at two evolution stages indicated by the core hy- 
drogen mass fraction (between curly brackets) . The unsta- 
ble radial orders (k m i n ,k max ) for different rotation rates 
f2 s (in units of f2 c ) and to- values are listed between round 
brackets . The numbers between square brackets list the 
oscillation periods (days) corresponding to k m i n and k max 
in the observer's frame. Positive periods correspond to 
prograde (in observer's frame) modes, negative periods 
to retrograde modes. Empty slots mean that no unstable 
modes could be detected. P s is the stellar rotation period. 
For X c = 0.4 and £l s ~ 0.4 no unstable modes were de- 
tected. 





Ps(d) 


m = 1 


m = 2 


m — 3 


m = 4 


0.4 


0.57 
{0.6} 




(15,18) 
[.396, .390] 


(15,19) 
[.224, .223] 


(17, 20) 
[.157, .157] 


0.3 


0.76 
{0.6} 


(15, 19) 
[3.6, 2.9] 


(14,21) 
[.54, .53] 


(14,21) 
[.301, .298] 


(17, 20) 
[.210, .210] 


0.2 


1.14 
{0.6} 


(13,22) 
[11., 6.] 


(14,22) 
[0.83,0.81] 


(15,20) 
[.453, .452] 




0.1 


2.28 
{0.6} 


(13,21) 
[79., 36.] 








0.3 


1.11 
{0.4} 






(25, 27) 
[0.438., 0.437] 


(25, 26) 
[.307, .307] 


0.2 


1.66 
{0.4} 


(21,27) 
[15.9, 10.9] 


(21,31) 
[1.22, 1.20] 


(21,30) 
[.662, .660] 




0.1 


3.33 
{0.4} 


(20,32) 
[121., 53.] 


(22,29) 
[2.48, 2.47] 







and the driving/damping zones, all as a function of tem- 
perature, for the slightly evolved 4M© MS model with 
X c = 0.60. The opacity Z-bump near 10 log T ~ 5.3 
coincides with the region where non-adiabicity of low- 
frequency, |cx/(2fi s )| < 1, oscillations becomes noticeable. 
A little further out the driving is, according to inequality 
H23fl . replaced by damping but the extra thermal energy 
exchange for relevant frequencies is small because of exist- 
ing strong non-adiabaticity. By comparing the dashed and 
dot-dashed curves it is seen that the thermal timescale Tth 
in the driving zone increases significantly when the star 
evolves further from the ZAMS and its surface layers be- 
come cooler. We determined the unstable n=0 g-modes 
with to=1 to to=4 in a 4M© main-sequence star rotat- 
ing at various angular speeds, for two evolution phases: 
X c =0.6 and X c =0.2, see Table Note that the (retro- 
grade) g-mode frequency |t7fc| increases with f2 s as the 
increasing Coriolis force renders the star 'stiffer' against 
oscillatory motions. In Fig. [S] it can be seen that, for a 
given rotation rate Q s , the n=0 g-mode frequencies \<Jk \ in- 
crease with increasing m- values, contrary to the behaviour 
of n = — 1 q-modes in Fig. |3| This is the usual g-mode 
behaviour, remember that n—0 means m=l in the non- 
rotating star terminology. The results in Table [3] show 
that in the little evolved MS model all combinations of 
Q s and to yield unstable g-mode intervals. However, for 
the evolved model with X c — 0.2 the opacity bump region 
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log_10(T) 

Fig. 4. MS star of 4M Q with X c =0.6 and Z=0.02. 
Full line: the modified Brunt- Vaisala frequency 
i^bv— sign(A/" 2 ) \J\N 2 1 , where TV is the Brunt- Vaisala 
frequency; dashed line: the 'adiabatic frequency' 
v a d — 2%/Ti ea k; dot-dashed line: again v a d, but now 
for X c = 0.2; dotted line: the value of the radial deriva- 
tive in expression (|23[l in arbitrary units, all versus 
temperature. 



has moved into the star with a corresponding longer local 
value for r t h- As a result the parameter combinations that 
correspond with higher mode frequencies, i.e. the larger m 
values and higher rotation rates, no longer yield unstable 
g-modes because the oscillations have become too adia- 
batic for destabilization by the K-mechanism. For f2 s =0.5 
all considered g-modes appear stable. 
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Fig. 5. Unstable n = retrograde g-modes: ak/(2S} s ) ver- 
sus radial node number k. For meaning of symbols see 
legend above, whereby 'mlo2' denotes case with m=l and 
fi s =0.2, etc. For given fl s the absolute value of the mode 
frequencies increase for increasing m(= /). Model: 
4M with X c =0.6 and Z=0.02. 



Table 5. Unstable n—Q retrograde g-modes for a 4 Mq 
star with Z= 0.02. See Table0]for explanation. 





P (d) 


TTl — 1 


m — 2 


771 —— 3 


Tfl = 4 


0.4 


0.63 
{0.6} 


(10,31) 

T 1 1 ro 1 

[-1., -158.J 


(10,31) 
[3.1, .64] 


(11,31) 
[.dz, .oo\ 


(12, 30) 

[ or ool 
[.oO, .22\ 


0.3 


0.84 
{0.6} 


(10, 30) 
[-1., -5.J 


(10,30) 

r a n i 1 1 
[-4.9, 1.1J 


(11,30) 
[1.7, .50] 


(12, 30) 

T 7fl QQl 

[. fU, .ooj 


0.2 


1.27 
{0.6} 


(10,29) 
hl.,-2.9] 


(11,29) 
[-1.7,27.] 


(11,30) 
[-2.4, .97] 


(12, 30) 
[-29., .62] 


0.1 


2.53 
{0.6} 


(10, 24) 
[-.99, -2.3] 


(10,27) 
[-.74,-5.1] 


(11,29) 
[-.68, -202.] 


(12, 30) 
[-.60,6.1] 


0.4 


1.45 
{0.2} 


(37, 44) 
[-3.3, -4.3] 








0.3 


1.93 
{0.2} 


(36, 53) 
[-3., -4.9] 


(37, 48) 
[11.5,4.5] 






0.1 


5.79 
{0.2} 


(30, 67) 
[-2.4, -5.6] 


(35,65) 
[-2.3,-8.8] 


(42, 49) 
[-2.8,-4.1] 





5.2.2. Unstable n = — 1 q-modes 

Let us now consider the behaviour of the n = — 1 q-modes 
in a 4Mq main-sequence star. We have extended the calcu- 
lations to a rather high angular rotation speed of tt s = 0.5 
for which the neglect of centrifugal distortion is question- 
able just to see the effect of stronger Coriolis forces on the 
oscillation. Fig. shows the eigenfunction X n (fi) for the 




Fig. 6. The angular eigenfunction Xn^d) versus /i = cos$ 
for four different m=2 modes, with respectively n = — 1 
and n=0, that bracket two instability domains for Q s = 
0.3, see Tables ISTol From top to bottom: Full line: q 12 \ long 
dashes: q2\\ short dashes: gio and dotted line: g3o (the 
latter shows the extra node and rotational confinement) 
all as listed in Tables [S] and The normalisation of y-axis 
is arbitrarily. Model: 4 M , X c = 0.6, Z=0.02. 

first and last unstable m=2 n=0 g-mode and m=2 n = — 1 
q-mode, bracketing their instability domains, as listed in 
Tables |S] and IH1 for fi s = 0.3. The two bracketing unstable 
g-modes are at the boundary, respectively well inside the 
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incrtial range: with 1/x ~ —1.09 and -0.60, respectively. 
For the latter low frequency g3o-mode rotational effects 
are thus significant and its eigenfunction X n exhibits the 
typical equatorial confinement of these modes: the eigen- 
function has significant amplitude only for \fi\ < 1/x. The 
two q-modes, by definition inside the inertial regime, have 
respectively the values 1/x ~ —0.302 and —0.271, but are 
seen to exhibit hardly any confinement toward low values 
of \i (these q-modes are odd and vanish at the equator). 
The q-modes, in spite of their weak compressions also un- 



Table 6. Unstable n = — 1 q-modes for a 4 Mq star with 
metal content Z— 0.02, see Table 01 for explanation. 





P s (d) 


m — 1 


m = 2 


m — 3 


m = 4 


0.5 


0.51 
{0.6} 


(10,27) 
[1.8,1.0] 


(11,25) 
[0.35,0.32] 


(12,22) 
[.198, .194] 


(13,21) 
[.139, .138] 


0.3 


0.84 
{0.6} 


(11,22) 
[4.9,2.5] 


(14,21) 
[0.60,0.58] 


(14, 17) 
[.333, .332] 




0.2 


1.27 
{0.6} 


(11,20) 
[13.4,6.1] 








0.1 


2.53 
{0.6} 










0.5 


0.74 
{0.4} 


(16,43) 
[2.6, 1.5] 


(16,41) 
[0.52,0.47] 


(17,39) 
[.291, .283] 


(19, 36) 
[.204, .202] 


0.3 


1.23 
{0.4} 


(15,39) 
[8.6,3.4] 


(16,35) 
[0.90, 0.84] 


(18,25) 
[.489, .486] 


(21,25) 
[.3416, .3412] 


0.1 


3.70 
{0.4} 


(19,23) 
[118,84] 








0.5 


1.16 
{0.2} 


(28,65) 
[4.1,2.4] 


(28,66) 
[0.81,0.74] 


(30, 65) 
[0.45, 0.44] 


(31,62) 
[.319, .316] 


0.3 


1.93 
{0.2} 


(27,66) 
[13.1,5.3] 


(28,63) 
[1.4, 1.3] 


(30, 58) 
[0.77, 0.75] 


(32, 53) 
[.534, .532] 


0.1 


5.79 
{0.2} 


(28,50) 
(253, 88] 









stable due to the K-mechanism, extent the instability re- 
gion down to lower frequencies than that of the normal g- 
modes of Table |SJ The q-modes have significantly smaller 
A-values and their angular eigenfunctions do not exhibit 
the strong rotational confinement to the stellar equator for 
small values of \l/x\ shown by normal retrograde g-modes. 

In Table we list the unstable q-modes with n = — 1 
and m=l to 4 at various evolutionary stages on the main- 
sequence and for angular rotation speeds up to fl s = 0.5. 
The r-mode frequencies are proportional to f2 s , and the 
same remains (roughly) true for the q-modes, although 
non-toroidal effects begin to become significant with in- 
creasing k. 

In order to understand the shifting boundaries of the 
instability regions with m and £l s one should again realize 
that for a given rotation speed f2 s the absolute value of 
the oscillation frequency \ak\ for n — g-modes increases 
with m(= I), see Fig. El while Fig.[7|shows that for n = — 1 
q-modes the reverse is true (for an explanation see section 
I5.1JI . It can be seen in Fig. e.g. the case with m = 1 
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Fig. 7. Unstable n = — 1 q-modes: <7fc/(2f2 s ) versus radial 
node number k. For meaning of symbols see legend above, 
whereby 'mlo2' denotes case with m=l and f2 s =0.2, etc. 
Contrary to the g-modes in Fig. [S] the absolute value of 
the mode frequencies |<7jt| decrease for increasing m for a 
given O s , see text. Model: 4M Q with X c =0.6 and Z=0.02. 

and f2 s =0.5, that sometimes not all radial orders in an 
instability domain are unstable according to our criteria 
(the mode amplitude should increase when the search for 
the resonant forcing frequency is made in the complex 
plane). In some cases the resonances remain too weak. 
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Fig. 8. Growth rates —Tm(a) versus A for unstable n 
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dotted squares, filled squares and dotted circles. Model: 
4Mq with X c =0.6 and Z=0.02 for Q s = 0.3. 

Fig. |H1 shows the imaginary parts of the complex forc- 
ing frequency a for which the resonant response in a AMq 
model is maximized for both the unstable n — — 1 q-modes 
and ?7=0 g-modes with m-values in the range 1-3. As ex- 
pected, for the quasi g-modes the highest growth rates 
are attained by the m=l modes which have the largest 
A values and highest oscillation frequencies a (see Fig. 
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and are thus closest in character to the g-modes. But in 
this little evolved stellar model (X c = 0.6) the driving 
layer is relatively close to the stellar surface and the (low- 
frequency) modes have rather small growth rates. 

In the more evolved stellar model with X c =0A and 
for the lowest rotation speed O s =0.1 the m=2 and to=3 
q-modes have oscillation frequencies too small for the k 
instability to be effective, the oscillations are too non- 
adiabatic in the K-bump region. In the most evolved stel- 
lar model (X c = 0.1), on the other hand, the sudden in- 
crease of Tad occurs further out in the star, see Fig.0] For 
ft s = 0.5 this results in stable m=l q-modes as the mode 
frequencies are in this case too high for efficient driving in 
the opacity bump region. 

Note again that for to > 2 the subsequent radial or- 
ders of the n = — 1 q-modes in an instability interval are 
densely spaced in period. For normal g-modes similar close 
spacing of radial orders k generally occurs for much larger 
values of k for which the oscillations are significantly more 
confined to the equatorial plane of the rotating star and 
for which the strong radiative damping prevents instabil- 
ity. 



(see Fig. HOf) . are depicted in Fig. [§| The eigenfunction 
Xn($) determines the latitudinal distribution of temper- 
ature, density (and small radial velocity) over the stel- 
lar surface. The shape varies little as a function of to. 

Table 7. Unstable n — —2 q-modes for a 3M Q MS star 
with Z=0.02, see Tabic 0] for explanation. For f2 s lower 
than listed no unstable modes could be found. 





P s (d) 


m = 1 


m = 2 


m — 3 


m = 4 


0.4 


0.57 
{0.6} 




(13,18) 
[0.32,0.31] 


(14, 19) 
[.206, .204] 


(16, 18) 
[.1509, .1507] 


0.3 


0.76 
{0.6} 




(14,17) 
[.434, .429] 


(16, 17) 
[.2764, .2760] 




0.4 


0.83 
{0.4} 


(19,25) 
[0.98, 0.95] 


(19,30) 
[0.47, 0.46] 


(21,31) 
[.301., .297] 


(22, 29) 
[.221, .220] 


0.3 


1.11 
{0.4} 


(20, 23) 
[1.36, 1.32] 


(20, 29) 
[0.64, 0.62] 


(21,28) 
[.405., .402] 


(24, 27) 
[.295, .295] 


0.2 


1.66 
{0.4} 


(21,22) 
[2.17,2.16] 


(20, 26) 
[.975, 0.964] 







5.3. Unstable n = —2 q-modes in 3-4 M & MS stars 

It appears that in the evolved 3-4M Q models even un- 
stable higher order n — — 2 quasi-toroidal modes are ex- 
cited by the K-mechanism. These modes have very small 
eigenvalues A_2 and very low oscillation frequencies, see 
Tables 17181 and Fig. ^| G-modes of similarly low oscilla- 
tion frequencies would be of much higher radial degree and 
have much larger A- values and be stable due to strong non- 
adiabaticity. In less evolved stellar models with X c > 0.4 
the n = —2 q-modes appear all stable. The angular eigen- 




0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1 
mu 

Fig. 9. The angular eigenfunction X n {$) versus /z for the 
unstable fc=17 n = — 2 modes at rotation rate Jl s =0.4: 
Solid line: m=2; dashed line: to=3 and dotted line: to=4. 
Model: 3 M Q , X c = 0.6, Z=0.02. 

function X n ($) for three (unstable) n = — 2 q 17 -modes 
with m—2, 3 and 4, all with frequencies \l/x\ < 0.115 



Table 8. Unstable n = —2 q-modes for a 4 M & star with 
X c — 0.2 and metal content Z= 0.02, see Table 0] for ex- 
planation. The oscillation periods in the corotating frame 
(between square brackets) are of order the rotation period 
over to due to the small values of |<7fc|. 





Ps(d) 


m = 1 


m = 2 


m = 3 


m = 4 


0.5 


1.16 
{0.2} 




(28, 43) 
[0.65,0.63] 


(29, 47) 
[0.418,0.410] 


(32, 47) 
[.306, .304] 


0.4 


1.45 
{0.2} 




(29, 40) 
[0.82, 0.80] 


(30, 44) 
[0.525,0.519] 


(34, 43) 
[.384, .382] 


0.3 


1.93 
{0.2} 




(32, 37) 
[1.106, 1.096] 


37 
[0.70] 





Figures and 0] show that in the evolved stars the 'adi- 
abatic' frequency rises steeply only near the outer edge 
of the opacity bump region. Low frequency q-mode oscil- 
lations apparently can still absorb sufficient heat during 
compression in the opacity bump zone, at least for rota- 
tion speeds f2 s > 0.3. When the opacity derivatives k t 
and k p are put to zero the instabilities disappear. 

For the higher stellar masses on the main sequence 
the opacity bump is closer to the surface and the non- 
adiabaticity too strong for destabilization of the very low 
frequency n = —2 q-modes. 

Fig. 1101 shows that the absolute value of the n — — 2 
mode frequencies (in corotating frame) \dk\ as a function 
of azimuthal index to behave differently from the n — — 1 
modes shown in Fig. [7| the absolute value of the frequen- 
cies of e.g. to=3 modes are higher than those of m=2, 
while for k > 16 the same holds for m=4 with respect to 
to=3. This can again be explained by considering equation 
(|2ip. For n — —2 the roots x" 1 have values —1/6, —1/6, 
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dial order k, for rotation rate fl s — 0.4: asterisks: m=4, 
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and Z=0.02. 
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Fig. 11. Growth rates — Zm(er) versus A for the unstable 
n = —2 q-modes shown in Fig.^] Symbols: m—2, 3 and 
4 are represented by +, x and *, respectively. Model: 3M Q 
with X c =0.6 and Z=0.02. 



—3/20 and —1/10 for m=l, 2, 3 and 4, respectively, i.e. 
lie closely together so that already for relatively small k 
we get the behaviour stated above. 

Fig. llll shows the imaginary parts of the complex forc- 
ing frequency a for which the resonant response is max- 
imized for the unstable n = —2 q-modes with m-values 
in the range 2-4. The real parts of the oscillation fre- 
quencies are shown in Fig. 1101 The very low frequency 
n = — 2 q-modes with small A- values have significantly 
smaller growth rates than the n = — 1 q-modes. The m=3 
modes attain the largest maximum growth rate. 



5.4. Unstable modes in a 6M MS star 

5.4.1. Unstable n=0 retrograde g-modes for Z=0.02 

Fig. 1121 shows the frequencies of unstable n=0 g-modes in 
the 6Mq model with X c =0.6, rotating at speeds O s =0.1 
and 0.3. By comparing with Fig.Elwhich depicts the same 
for a AMq stellar model it appears that in this more mas- 
sive star the instability intervals are less extended and 
terminate at larger fc-values (Afc ~ —10), while the high 
frequency edge is shifted to slightly lower k with Ak ~ — 2. 
This is caused by the shorter thermal timescale in the driv- 
ing zone of the more massive star. Apart from the slight 
shift in fc, the frequencies of the unstable g-modes are 
very similar for the 4M and 6M Q models. Table [§] lists 
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Fig. 12. Unstable n = retrograde g-modes: cjfc/(2f2 s ) 
versus radial node number k. For meaning of symbols, 
see Figure, whereby 'mlo2' denotes case with m=l and 
f2 s =0.2, etc. For given J7 S the absolute value of the mode 
frequencies \ak\ increase for increasing m=l. Model: 6M Q 
with X r =0.6 and Z=0.02. 



the g-mode instability intervals for various rotation rates 
and evolution stages of the 6M© MS star. The table shows 
that we expect to observe many unstable g-modes with 
a retrograde character in observer's frame (negative peri- 
ods). By writing a = a + m fi s we can understand that the 
modes become prograde in the observer's frame for larger 
m-values and or higher stellar rotation rates as listed in 
table IH1 The unstable g-modes, both prograde and ret- 
rograde, occur in many wide instability intervals which 
is, however, not observed in SPB stars. We list here only 
the n=0 g-modes, but there exist also many higher order 
(n > 0) unstable g-modes. It is a serious puzzle why all 
these apparently unstable g-modes are not observed. 

5.4.2. Unstable n = -1 q-modes for Z=0.02 

The frequencies of unstable n = — 1 q-modes detected 
for the 6M Q model with X c =0.2 are shown in Fig. 



The instability regions (k 
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Table 9. Unstable n—0 g-modes for a 6Mq MS star with 
Z=0.02, see Table 0] for explanation. 





P (d) 


771 — 1 


m = 2 


vn — 3 


0.3 


0.98 
t u - D f 


(9, 18) 

r i ^ 9 si 

[ 1.6, — Z.SJ 


(9, 18) 

f 1 R 1 71 


(9,20) 

19 fi^l 


0.2 


1.47 
10 61 


(10, 14) 
[— 1 4 —1 91 


(9, 16) 
1— 1 9 28 1 


(9, 18) 
r_9 r i 7i 


0.1 


2.53 
10 61 


(10, 24) 
[_ qq _o 31 


(10,27) 
1 — 74 —5 11 


(11,29) 
[_ r« —202 1 


0.5 


0.85 
{0.4} 


(12,37) 
r i e ii ni 

[—1.0, 11. UJ 


(12,37) 

fn c n Ql 
[Z.O, U.oJ 


(13,38) 

fn 7 n /ii 
[UJ, U.4J 


0.3 


1.55 
{0.4} 


(13,31) 
[-1.7,-5.0] 


(13,32) 
[-6.6,2.1] 


(13, 34) 
[5.5, 0.89] 


0.1 


4.26 
{0.4} 


(14,20) 
[-1.8,-2.5] 


(13, 25) 
[-1.2,-3.5] 


(13,29) 
[-0.95, -6.1] 


0.5 


1.32 
{0.2} 


(21,68) 
[-2.5, 13.2] 


(21,60) 
[4.0, 1.2] 


(21,59) 
[1.2,0.67] 


0.3 


2.20 
{0.2} 


(21,61) 
[-2.4,-10.4] 


(21,62) 
[-6.5,2.9] 


(21,63) 
[31.0, 1.3] 


0.1 


6.59 
{0.2} 


(21,42) 
[-2.3, -4.7] 


(21,50) 
[-1.7,-7.6] 


(21,57) 
[-1.3,-21.1] 




log_T0(T) 

Fig. 13. MS star of 6M Q with X c =0.6 and Z=0.02. 
Fw/Z line: the modified Brunt- Valsala frequency 
vbv= sign(A/" 2 ) a/| A/" 2 1 , where TV is the Brunt- Valsala 
frequency; dashed line: the 'adiabatic frequency' 
z^ ac i = 2 7r/r; ea fe; dot-dashed line: again ^ a( j, but now 
for X c = 0.4; dotted line: the value of the radial deriva- 
tive in expression (|23|l in arbitrary units, all versus 
temperature. 



the g-modes: due to the lower oscillation frequency of q- 
modes the low-frequency (high k) end is less extended as 
the opacity bump region becomes too non-adiabatic. Note 
that all unstable q-modes in Fig. 1141 have \l/x\ < 0.45, 
while many of the unstable g-modes in the same stellar 
model (except the m=l modes) lie outside the inertial 
range, see Fig.1121 In Table ITU1 we list the unstable n = — 1 
modes with m—1, 2 or 3 for a 6 Mq main-sequence star. 
The frequency of the m=4 q-modes appears too low for 
destabilization by the n mechanism. With increasing stel- 



lar mass the opacity bump driving region for the insta- 
bilities moves closer to the stellar surface and the steep 
rise of T a d occurs in the 6Mq model at higher tempera- 
ture, see Fig. El For the unevolved 6 Mq models (with 
X c > 0.4) the thermal timescale in most of the opacity 
bump region is therefore shorter than the oscillation pe- 
riods (in corotating frame) of the q-modes, except for the 
fastest rotation rates considered. For this reason we find 
only m—1 instabilities and those only for f2 s > 0.4. When 
the star evolves away from the ZAMS the effective tem- 
perature decreases and the driving zone moves inwards, so 
that the local thermal timescale increases. For the models 
with X c < 0.4 we therefore find more unstable q-modes 
and also instabilities at lower rotation rates, although not 
for f2 s < 0.3 as in a 4 Mq star: the frequencies are in that 
case too low. 
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Table 10. Unstable n = — 1 q-modes for a 6M Q MS star 
with Z=0.02, see Table 0] for explanation. No unstable 
m=4 modes were found. 





Ps(d) 


m = 1 


m = 2 


m = 3 


0.5 


0.85 
{0.4} 


(15,23) 
[2.7,2.1] 






0.4 


1.07 
{0.4} 


(16,20) 
[4.0, 3.3] 






0.5 


1.32 
{0.2} 


(22,49) 
[4.9, 2.9] 


(24, 43) 
[.92, .87] 


(27,37) 
[0.516,0.511] 


0.4 


1.65 
{0.2} 


(22,44) 
[7.8, 4.4] 


(25, 38) 
[1.16, 1.12] 




0.3 


2.20 
{0.2} 


(24, 39) 
[13.6, 8.0] 
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5.4.3. Unstable n = — 1 q-modes for Z=0.03 

When the metal fraction is increased to Z=0.03 the opac- 
ity bump effect gets stronger and more q-modes become 
unstable, as can be seen by comparing Table ^| with 
Table ^2 With the enhanced opacity bump the higher 
m modes can also be destabilized, for high rotation rates 
ft s >0.3 even m=4 in the evolved X c —0.2 stellar model. 

Table 11. Unstable n = — 1 q-modes for a 6Af© MS star 
with Z=0.03, see Table 0] for explanation. 





P s (d) 


m = 1 


m = 2 


m = 3 


m = 4 


0.5 


0.93 
{0.4} 


(14, 32) 
[3.4, 2.0] 


(15,27) 
[.65, .61] 


(17, 24) 
[.365, .361] 




0.4 


1.17 
{0.4} 


(14, 28) 
[5.3,2.9] 


(16,25) 
[.82, .79] 






0.3 


1.55 
{0.4} 


(15,26) 
[9.3, 5.3] 








0.5 


1.47 
{0.2} 


(22, 56) 
[6.0,3.1] 


(24, 58) 
[1.03,0.94] 


(25, 53) 
[.578, .564] 


(27, 47) 
[.406, .403] 


0.4 


1.84 
{0.2} 


(22, 60) 
[9.6,4.2] 


(24,54) 
[1.31, 1.22] 


(26,47) 
[.726, .714] 


(32, 37) 
[.508, .507] 


0.3 


2.45 
{0.2} 


(23, 53) 
[18., 7.2] 


(25,46) 
[1.77, 1.69] 


(28, 34) 
[.971, .968] 





Fig. 1 151 shows the imaginary parts of the complex forc- 
ing frequency a for which the resonant response is maxi- 
mized for both a set of n = — 1 q-modes and n=0 g-modes 
for various m- values. Again the highest frequency m=l q- 
modes show the largest growth rates. It can be seen that 
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Fig. 15. Growth rates — Xm(cr) versus A of unstable n = 
— 1 q-modes and n=0 g-modes. Symbols: q-modes with 
m=l,2 and 3: +, x and *; g-modes with m=l,2 and 3: 
dotted squares, filled squares and dotted circles. Model: 
6M Q with X c =0.2 and Z=0.03 for Q s = 0.3. 



5.5. Unstable n = — 1 modes in a 8M MS star with 
Z=0.03 -> LPV in Be stars? 

Fig. 1161 shows the driving- and damping regions in a 8 
M Q MS star with X c = 0.2 and Z=0.03. It can be seen 
that for this more massive star the opacity bump region is 
convectively unstable. Nevertheless, by assuming the con- 
vective flux is not heavily perturbed by the oscillation, we 
have sought for unstable modes, see Table IT2l We have ex- 
tended the calculations up to completely illegal rotation 
speeds for which the (neglected) rotational deformation of 
the star becomes strong in order to compare our results 
with observations of rapidly rotating Be-stars. It is tempt- 




7.5 7 6.5 6 5.5 5 4.5 



log_10(T) 

Fig. 16. MS star of 8M with X c =0.2 and Z=0.03. 
Full line: the modified Brunt- Vaisala, frequency 
^bv= sign(A/' 2 )-\/|A/' 2 |, where J\f is tfhe Brunt- Vaisala 
frequency; dashed line: the 'adiabatic frequency' 
v a( i — 2-K/rieak] dotted line: the value of the radial 
derivative in expression (|23|l in arbitrary units, versus 
temperature. 



Table 12. Unstable n = -1 q-modes in a 8 M© MS star 
with Z = 0.03, see Tabled for explanation. 





Ps(d) 


m = 1 


m = 2 


m = 3 


0.9 


0.90 
{0.2} 


(20,34) 
[2.21,1.70] 


(21,44) 
[.592, .552] 


(23,39) 
[.345, .334] 


0.7 


1.15 
{0.2} 


(21,37) 
[3.3, 2.4] 


(22,38) 
[.78, .74] 


(25,33) 
[.446, .442] 


0.6 


1.34 
{0.2} 


(21,36) 
[4.3, 3.0] 


(23,35) 
[.92, .88] 




0.5 


1.61 
{0.2} 


(21,34) 
[6.0,4.1] 


(25,31) 
[1.11, 1.09] 




0.4 


2.02 
{0.2} 


(23, 30) 
[8.7, 6.6] 







the maximum growth rate of m=l q-modes becomes al- 
most comparable (<~ 20%) to that of n—0 g-modes in this 
more evolved (X c = 0.2) main-sequence model. 



ing to identify the observed line profile variability (LPV) 
in the well observed Be star /j, Cen (Rivinius et al. |2001 ) 
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in terms of n = — 1 q-modc oscillations in a rapidly rotat- 
ing B-star. The observed variations are indeed all prograde 
in the observer's frame and fall in two narrow period in- 
tervals of about 0.50 days (m=2) and 0.28 days (m=3). 
Qualitatively this is exactly what is found for the unsta- 
ble n =-1 q-modes calculated here (see Table H"2*|l . although 
the exact periods do not match. But that is not surpris- 
ing in view of the over-simplistic modelling of extremely 
rapidly (differentially) rotating Be-stars. 

6. Summary and discussion 

In this paper we study the stability of quasi g-modes ('q- 
modes'), a branch of oscillation modes that exists in be- 
tween the normal retrograde g-modes and the r-modes in 
rotating stars. It appears that these retrograde (in corotat- 
ing frame) low frequency oscillations can be destabilized 
by the K-mechanism in the driving region connected with 
the opacity bump close to the stellar surface. From rela- 
tively low angular rotation speeds upwards the unstable 
q-modes with m > 1 occur in a few narrow period bands 
(for different m- values) in which the different radial orders 
lie densely packed together. For m=l the periods of q- 
modes are typically larger than a day, up to many days for 
low radial orders. In the observer's frame all q-modes are 
prograde. The predicted observable q-mode spectrum is 
thus in sharp contrast to the predicted spectrum of unsta- 
ble g-modes which occur in much wider period intervals, 
whereby modes of different radial order k should be sep- 
arately observable with even modest time resolution. The 
unstable q-modes seem to fit the observed rather scarce 
frequency spectrum in SPB stars and [3 Cephei stars (e.g. 
Stankov ct al. 2002) better than normal g-modes. The ob- 
served closely spaced frequencies in e.g. the SPB star HD 
160124 (Waelkens, 11991(1 cannot be explained by unsta- 
ble g-modes, but is typical for the q-mode spectrum. It 
is, however, doubtful whether the rotational confinement 
of g-modes towards the stellar equator in rotating stars 
can explain away the many predicted unstable g-modes 
(see Townsend IMIS) . It remains mysterious why we don't 
observe the many predicted unstable g-modes and would 
instead observe the (weaker) unstable q-modes. It helps 
in any case that the unstable q-modes hardly have a ten- 
dency to focus on the equatorial region of rotating stars. 

Although the more massive j3 Cephei pulsators seem 
to show clustered periods similar to SPB stars (but with 
smaller periods) they cannot be explained straightfor- 
wardly in terms of q-modes because in stars more massive 
than 8Mq the thermal timescale in the the Z-bump driv- 
ing region becomes too short to enable destabilization of 
the q-modes. Q-mode instability for these more massive (3 
Cephei pulsators would require the Z-bump region to be 
deeper inside the star than predicted by the current stellar 
evolution models. 

Nevertheless, the observed LPV in Be stars like fj, Cen 
is (qualitatively) consistent with the calculated properties 
of q-modes in a rapidly rotating 8M0 B-star: two narrow 
period intervals corresponding to two unstable series (for 



m=2 and m=3) of closely spaced q-modes, see previous 
section. One may speculate about the cause of the related 
outbursts in this Be-star in terms of q-modes. Excitation 
of the (negative angular momentum) retrograde q-modes 
in the star will accelerate the rotation of the stellar mate- 
rial in the driving region just beneath the stellar surface, 
perhaps bringing the surface layers up to ejection veloci- 
ties. A good understanding of the Be phenomenon requires 
a much more detailed study with differential rotation in- 
cluding the centrifugal force and non-linear effects. 
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